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ABSTRACT 

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (M* > 10 11 Mq) 
galaxies at 2 < z < 3 in the GOODS (Great Observatories Origins Deep Survey) North field using 
the Spectral and Photometric Imaging Receiver (SPIRE, Griffin et al. 2010) instrument aboard 
the Herschel Space Observatory. To conduct this analysis we take advantage of the data from 
the HerMES key program. The sample comprises 45 massive galaxies with structural parameters 
characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with 
Spitzer 24fim data, as a function of the morphological type, mass and size. We find that 26/45 sources 
are detected at MIPS-24^m and 15/45 (all MIPS-24/jm detections) are detected at SPIRE-250^m, 
with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star 
formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into 
account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that 
the mean SFR for the spheroidal galaxies (~ 50 — lOOMQyr -1 ) is substantially (a factor ~ 3) lower 
than the mean value presented by disk-like galaxies (~ 250 — SOOMoyr^ 1 ). 



Key words: infrared: galaxies - galaxies: evolution - galaxies: high-redshift - galaxies: star for- 
mation 
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1 INTRODUCTION 

One of the most intriguing recent results in extragalactic 
astrophysics is the discovery that massive galaxies (M > 
10 n M Q ) at high redshift were on average more compact 
than their local counterparts (Daddi et al. 2005; Trujillo 
et al. 2006, 2007; Longhetti et al. 2007, Buitrago et al. 
2008, B08 hereafter). Although the formation mechanism 
for such compact galaxies is not yet perfectly established 
(e.g. Khochfar & Silk 2006, Ricciardelli et al. 2010), the 
proposed scenario that drives these compact galaxies into 
the local population of massive ellipticals likely includes dry 
major mergers or, the currently favoured scenario, minor 
mergers with less dense galaxies (Naab et al. 2007; Hopkins 
et al. 2009). 

In order to shed light on how these compact galaxies 
have formed the bulk of their mass, the characterisation of 
their star formation rate (SFR) history is required. While 
their morphological properties seem to be quite simple, their 
intrinsic physical properties are still not well characterised. 
Different groups have found contrasting results on the corre- 
lation between the SFR and structural properties. Kriek et 
al. (2009), for example, claim a striking correspondence be- 
tween the SFR and structural parameters (sizes, morpholo- 
gies), with the compact red galaxies forming stars at a rate 
one order of magnitude lower than larger blue galaxies. In 
contrast, Perez-Gonzalez et al. (2008, PG08 hereafter) using 
extremely deep 24 /im imaging of the Extended Groth Strip 
found that at z ~ 2, ~ 80% of the disk-like galaxies and 
~ 50% of the compact spheroid-like galaxies are detected, 
with similar mean SFRs of 6O-lOOM yr _1 . However, within 
each morphological category the size is not related to the 
24 fim detection. 

In the present work we focus on a sample of massive 
(M* ~ 10 11 - 1O 12 M ) high redshift (z ~ 2 - 3) galaxies 
having deep HST data with sizes and structural properties 
measured by B08 and examine their FIR properties, as mea- 
sured from _ffersc/ie/Q(Pilbratt et al. 2010) within the Her- 
schel Multi-tiered Extragalactic Survey (HerMES, Oliver et 
al. 2010b), as a function of morphology and compactness. 
This paper represents the first exploratory step in this di- 
rection. 

Throughout this paper we assume a cosmology with 
H = 70kms- 1 Mpc-\ Q M = 0.3, fi A = 0.7. 



2 THE DATA 

The parent sample of galaxies originates from the work of 
B08. It covers the GOODS North and South fields and has 
been imaged as part of the GOODS NICMOS Survey (GNS; 
PI C. Conselice). We make use only of the Northern sub- 
sample in this work. 

We take advantage of the HerMES science demonstra- 
tion phase (SDP) observations carried out at 250, 350 and 
500 /im, which cover a large (~ 30 arcmin by 30 arcmin) area 
centred in the GOODS-N field. The instrument characteris- 
tics and capabilities are described in Griffin et al. (2010). 

1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 



2.1 The GNS sample 

The GNS is a large HST NICMOS-3 camera program of 60 
pointings centred around massive galaxies at z ~ 1.7 — 3 in 
the F160W (H) band in the GOODS North and South fields. 
The pointings were optimized to observe as many high-mass 
M* > 10 11 Mq galaxies as possible, with the selection of 
these targets described in Conselice et al. (2010). 

Stellar masses and photometric redshift estimates 
for this sample were derived exploiting multi wavelength 
(BVRIizJHK) data from GOODS (Giavalisco et al. 2004), 
using Bruzual & Chariot (2003) stellar population synthesis 
models, assuming a Chabrier (2003) Initial Mass Function 
(IMF). Only four spectroscopic redshifts are available for 
GOODS-N (Barger et al. 2008), showing a mean Az/(l + 
z) ~ 0.08 with respect to the photometric redshifts. Details 
of SED fitting techniques are in Conselice et al. (2007). As 
explained in the following section, in deriving the SFR we 
are implicitly assuming a Salpeter (1955) IMF. Therefore, 
for consistency and in order to derive the specific star forma- 
tion rates (sSFR), we report all the masses used in this work 
to a Salpeter IMF. To scale the stellar masses of the par- 
ent catalogue from Chabrier to Salpeter we used the follow- 
ing equation: logA/*(Salpeter)=logM*(Chabrier)+0.23 and, 
when needed (for the Shen et al. 2003 stellar mass-size 
relation), we also scale from the Kroupa (2001) IMF to 
Salpeter using: logAf*(Salpeter)=logA4*(Kroupa)+0.19 (see 
e.g. Cimatti et al. 2008). Structural parameters and sizes 
for the GNS sample were measured by using the GALFIT 
code. The light distribution has been modelled with a Sersic 
model, deriving Sersic indicx n, axis ratio b/a and effective 
radius along the semi-major axis a e . Our measured sizes are 
circularized so that ro = a t \/l — e with e the projected el- 
lipticity of the galaxy. Details on the fitting procedure and 
reliability of the measurements can be found in B08 and 
Trujillo et al. (2007). 

This data set, containing 45 galaxies for GOODS-N, 
represents the largest sample of high-z massive galaxies to 
date, that have high-quality HST (NICMOS) data and mea- 
sured structural parameters. This allows us to look for the 
first time at the possible correlations between rest-frame far 
IR properties, exploiting the available Herschel data, and 
the structural parameters (sizes and Sersic indices) of these 
high-z massive galaxies. We will refer to this catalogue as 
the parent catalogue throughout the paper. 

2.2 The Herschel data 

The SDP GOODS-N observations are among the deep- 
est possible with SPIRE, since the instrumental noise is 
lower then the confusion noise from overlapping background 
sources. In building the catalogues two different approaches 
have been used. In one approach blind source extraction 
was performed, resulting in single-band catalogues at 250, 
350 and 500/im (Smith et al. 2010, in prep.). In the second 
approach, flux densities and S/N ratios were obtained from 
PSF fitting on 24pm priors using a new source extraction 
method developed for HerMES called 'XID' (for details see 
Roseboom et al. 2010). Cross identifications are performed 
in map space so as to minimize source blending effects and 
allow the recovery of the faintest sources, close to the con- 
fusion limit. 
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Figure 1. Stellar mass-size relation for disk-like (left panel) and spheroid-like (right panel) galaxies. The dashed line in both panels 
indicates the local relation (Shen et al. 2003). The shaded area represents the stellar mass-size relation for galaxies between z = 2 (upper 
envelope) and z = 3 (lower envelope) from B08. Black circles represents galaxies in the parent catalogue, all filled circles are 24/im 
detections, while coloured filled-circles (blue and red for disk-like and spheroid-like galaxies respectively) indicate the 250/im detection. 



Within a field of approximately 10 arcmin by 15 arcmin, 
which has been observed with Spitzer and HST, there are 
a total of ~1500 24 fj,m sources which meet the reliability 
criterion of S/N^5 and 524 ^ 20/iJy. Herschel detects about 
a third of them in at least one band (as reported by Elbaz 
et al. 2010). 

Following the approach of Elbaz et al. (2010), we use 
SPIRE measurements down to the 5a limits of the prior 
catalogue of 4.4 mjy at 250/an. We note that these mea- 
surements lie close to the SPIRE confusion noise of 5.8 mjy 
(Nguyen et al. 2010). However, this limit is a spatially av- 
eraged statistical limit which considers that galaxies are ho- 
mogeneously distributed in the field and all affected in the 
same way by close neighbours. We take advantage of the 
higher resolution 24 /jm image to flag 'clean' galaxies (for 
which SPIRE flux densities can be potentially more robust) 
by requiring that that there is at most one bright neigh- 
bour within 20 arcsec (close to the full width half maxi- 
mum, FWHM, of the SPIRE 250 /im band) with 24 /im flux 
higher than 50 per cent of the central 24 fj,m counterpart. 
By matching the resulting catalogue with the GNS cata- 
logue from B08 we obtain a sample of 15 detected sources in 
total. Additionally we have 11 24/im priors without a 'clean' 
SPIRE detection. Of the 15 detections, 10 galaxies show de- 
tection at 350 /im and only 3 have detection at 500 /jm. Due 
to the large PSF of SPIRE and the associated deblending 
problems, especially at longer 350 and 500 /im wavelengths, 
we focus our attention in this work on the 250 /jm detections. 
As an additional check we also looked for possible 250/mi de- 
tections not associated with the prior 24/mi sources. We can 
confirm that none of the 19 GNS galaxies without 24 /im 
detections have a clear detection at 250 /im. 



3 ANALYSIS AND RESULTS 



3.1 IR luminosities and SFR 

We have performed an SED fitting procedure in order to de- 
rive the total IR luminosity of each source, by comparing the 
optical-to-IR photometry with a library of template SEDs 
of local objects from Polletta et al. (2007) and adding a few 
modified templates (Gruppioni et al. 2010). We have sep- 
arately analyzed objects with SPIRE and/or MIPS detec- 
tions. When limited to the 24 /im band, large extrapolations 
are required. However, with our technique (see Rodighiero et 
al. 2010) the inclusion of the whole SED in the fitting proce- 
dure allows us to fully exploit the photometric information 
and to determine appropriate k-corrections. As described in 
Rodighiero et al. (2010), the advent of Z/erschel data con- 
strains the far-IR normalization of the SED, providing more 
accurate bolometric luminosities (and consequently SFR val- 
ues). In both cases, we forced the spectral fit to reproduce 
the far-IR Spitzer and Herschel data points by weighting 
more the FIR points. 

As a measure of the total IR luminosity, our adopted 
procedure has been to integrate for each object the best-fit 
SED, in the rest-frame [8-1000/im] interval. For the galaxies 
detected with SPIRE or at 24 /im, the instantaneous SFR 
was then estimated using the calibration of Kennicutt et al. 
(1998): SFRIMoyr" 1 ] = 1.7 x 10" 10 L([8 - lOOO/xm])/L , 
under the assumption that the bulk of the bolometric emis- 
sion is dominated by star-formation with a Salpeter IMF. 
We have also checked that the galaxies in our sample do not 
have X-ray counterparts, ensuring the absence of optically 
thin AGNs detected by the Chandra Space Observatory. 
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Table 1. SPIRE 250 fim and MIPS 24 /im detection fraction, mean SFR and sSFR. 









SPIRE 250 /im 






MIPS 24 nm 






Ntot 


Detection 


SFRIMoyr- 1 ] 


sSFR[Gyr _1 ] 


Detection 


SFR^Moyr- 1 ] 


sSFRIGyr- 1 ] 


All 


45 


33 ± 8 % 


240 ± 50 


1.05 ± 0.33 


58 ± 11 % 


225±60 


0.95 ± 0.25 


n> 2 


16 


13 ± 9 % 


95 ± 35 


0.34 ± 0.13 


38 ± 15 % 


45 ±25 


0.17 ± 0.11 


n< 2 


29 


45 ± 12 % 


310 ± 45 


1.25 ± 0.25 


69 ± 15% 


280±70 


1.22 ± 0.30 



3.2 Results 

We split our sample into early- (or spheroid) and late-type 
(or disk-like) galaxies, according to their measured Sersic 
indices in H band from the NICMOS images. B08 showed 
that a value of n — 2 is optimal to separate galaxies at these 
high redshifts into the two morphological classes when us- 
ing NICMOS data. Therefore, we define spheroids as galax- 
ies with n > 2 and disk-like galaxies those with n ^ 2. We 
obtain for the parent catalogue 29 disk-like galaxies and 16 
spheroids (that is ~ 65 per cent of the sample is composed 
of disk-like galaxies, while ~ 35 per cent are spheroids). 
We have also checked, through one and two dimensional 
Kolmogorov-Smirnov tests, that the underlying distribution 
of stellar masses and redshift are compatible with being ex- 
tracted from the same parent distribution (at a 99 per cent 
of confidence level). This should insure the detection frac- 
tions are not largely biased by the underlying distributions. 

As an additional criterion, we split our sample according 
to the surface mass density (defined as D — 0.5 ■ M/{-ktq)) 
of our galaxies, adopting as threshold the value of D t = 
9 x lO 9 M0kpc -2 , that is the median value for our sample. 
The resulting fraction of MIPS detected over-dense (with 
respect to the median value) galaxies result to be ~ 41±13%, 
while for under-dense galaxies is ~ 74 ± 18%. Analogously, 
the SPIRE detected fractions result to be ~ 18 ± 9% and 
~ 48 ± 14% for over-/under-dense galaxies respectively. In 
both cases, relatively less dense galaxies result to have higher 
detection fractions. 

We have also checked through the Fisher one-sided ex- 
act test that the detected fractions are in agreement with 
the hypothesis that there is a connection between detection 
and morphology (at a ~ 98 per cent confidence level). 

As a first result in Fig[T] we show the stellar mass-size 
relation for disk-like (left panel) and spheroid (right panel) 
galaxies, including MIPS and SPIRE detections. The dashed 
lines in each panel indicate the local mass-size relation (Shen 
et al. 2003) . The shaded area represents the region enclosed 
between the mass-size relation for high redshift galaxies at 
2 iS z iS 3 following the evolution found in B08. All circles 
represent galaxies in the parent catalogue, with filled circles 
(black and coloured) being 24/^m detections and those in 
colour (blue and red for disk-like and spheroid-like galaxies 
respectively) being the 250pim detections. Both relations, lo- 
cal and high redshift, have been transformed to the Salpeter 
IMF using the equations reported in Section 2.1. 

We note that SPIRE detections correspond mainly to 
larger and more massive disk-like galaxies (~ 45 per cent 
of the sample has a 250^tm counterpart, and ~ 70 per 
cent has a MIPS detection, see Fig[T}, while spheroids 
are mostly undetected (only ~ 13 per cent have 250/im 
detection and only ~ 38 per cent are MIPS detected). 



The mean radius for 250/im-detected disk-like galaxies in 
our sample is log(ro)=0.50 ± 0.10, while for non-detected 
sources it is log(ro)=0.37 ± 0.08. Furthermore, the mean 
mass for 250/im-detected disk-like galaxies in our sample 
is log(Mass) = 11.45 ± 0.10, while for non-detected sources it 
is log(Mass) = 11.40 ± 0.05. This results reflect in the trends 
shown in Fig. 2 and are in agreement with previous sub-mm 
studies where it is found that the radius of the 850/im- 
detected sources is larger than average optically-detected 
galaxies at essentially all redshifts (e.g. Chapman et al. 2003, 
Pope et al. 2005). In the case of spheroid- like galaxies a 
higher detection fraction is also suggested for larger objects, 
but this is not so clear, due to the small statistics. The 
two spheroids detected with SPIRE have a Sersic index just 
above the threshold value (2.2 and 2.35). These are classified 
using a fit to the NICMOS imaging, but from inspecting the 
ACS images there is at least some hint of possible structure 
associated with one of the two spheroids. Additionally, we 
note that all the galaxies with a mass-size relation in bet- 
ter agreement with the local one (see also Fig[T] left panel) 
could correspond to interacting systems, as suggested by the 
inspection of ACS images. If confirmed, large sizes would 
not indicate a big galaxy, but rather a compact galaxy in 
a merging phase. It is beyond the goals of the present Let- 
ter to examine the details of these objects, that are instead 
deferred to a future work. 

In the four panels of Fig[2] we show the trends for the 
detection fractions for MIPS- and SPIRE-detected galaxies 
as a function of Sersic index (upper left panel), effective ra- 
dius (upper right panel, units of kpc) , stellar mass (bottom 
left panel, in Log(M0)) and redshift (bottom right). The 
points represent the centres of three equally populated bins 
in the parent sample. The black dots indicate the mean de- 
tection fractions with Poissonian error bars for MIPS, while 
the red dots are the fractions of SPIRE 250/im detected 
sources. These panels indicate a trend in the sense of higher 
detection for large, massive, disk-like galaxies. 

In Table [T] we summarise the results for the detection 
fraction, SFR and sSFR by dividing the sample according 
to morphology, for MIPS- and SPIRE-detected galaxies. The 
fractions of MlPS-detected galaxies are compatible with pre- 
vious results (e.g. PG08). 

The mean SFR and sSFR values are derived by taking 
into account non-detections for SPIRE, as explained below. 
We assume as lower limit null SFRs for 250/^m non-detected 
galaxies. Upper limits are derived by assuming a flux in the 
250^tm band equal to the 5cr limit of 4.4mJy for all non- 
detected sources and performing a fit of the mean SED for 
the whole sample (and the disk-/spheroid-like galaxies sub- 
samples) of non-detections. The reported errors refer to the 
semi-interval between the upper and the lower limits. 
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When computing the MIPS SFR/sSFR values we take 
into account the systematic effects measured in recent Her- 
schel works (see Elbaz et al. 2010; Nordon et al. 2010; 
Rodighiero et al. 2010), where it is found that MIPS based 
LIR estimates suffer of a systematic overestimation effect. 
The mean value for the ratio between 24pm and 250/im 
based LIR measurements for our sample results to be ~ 3, 
in agreement with the value found by Elbaz et al. (2010; see 
their Fig. 2) for high-z galaxies and using PACS+SPIRE 
data. Nordon et al. (2010) report a higher value (from 4 to 
7 for high-z galaxies) but this could be related with the fact 
that in this latter work the factor is derived by fitting purely 
PACS data. 

The mean values for the SFR of the sub-samples defined 
according to the morphological type indicate higher SFR 
for disk-like galaxies with respect to spheroidal-like galax- 
ies. This behaviour is reflected in the mean sSFR values de- 
termined for the two sub-samples. Furthermore, the SFR 
estimates well match other literature results (e.g. PG08, 
Rodighiero et al. 2010, Viero et al. 2010), once we restrict 
to the same mass-redshift range and convert the results to 
the same IMF. In a forthcoming Paper we will perform a 
complementary full stacking analysis (and take advantage 
of the larger statistics provided by the full GNS sample) in 
order to derive more robust estimates. 

In Fig(3] we show the mean specific SFR of massive 
galaxies detected by Herschel, as a function of redshift and 
morphology, compared to PG08, Oliver et al. (2010a) and 
Rodighiero et al. (2010). When segregating the sample based 
on the Sersic indices, PG08 find that the specific SFRs of 
massive (M > 10 11 Mq) spheroid galaxies detected at 24/im 
evolve as (1+jz) 5,5±0,6 from z = to z = 2 (red dashed line in 
FigEJi while the evolution for disk-like galaxies detected at 
24^m goes as (1 + z ) 3 - 6±0 - 3 (blue dotted line in Fig[3]). Also 
showed are the mean relations for early- (dark green dashed 
line) and late-type (light green) galaxies obtained by Oliver 
et al. (2010a) analysing data from the Spitzer Wide-area In- 
fraRed Extragalactic Legacy Survey (SWIRE; Lonsdale et 
al. 2003). The red dot indicates the mean sSFR for spheroid- 
like galaxies detected with SPIRE, while the blue dot is for 
the disks-like galaxies. The black point is the global mean 
value (the points are slightly shifted in z for clarity). The yel- 
low shaded region indicates the trend found by Rodighiero 
et al. (2010) using only Herschel-PACS (Photodetector Ar- 
ray Camera and Spectrometer, Poglitsch et al. 2010) data 
for a larger sample of GOODS-N galaxies (in the mass range 
10 n M o < M < lO n ' 5 M ). 



4 DISCUSSION AND CONCLUSIONS 

In this paper we have analyzed the FIR properties of an 
optically selected sample of 45 high redshift (1.7 ^ z ^ 3) 
massive (M > 10 11 Mq) galaxies as a function of their struc- 
tural parameters. Our analysis is based on the measurement 
of the specific SFR for each galaxy based on their IR emis- 
sion, exploiting the deep _ffersc/ie/-SPIRE data obtained by 
the HerMES Project in the GOODS North field. Our main 
conclusions are the following: 

• we globally detect 33 per cent of the galaxies of the 
parent GNS sample at least in one SPIRE band, this is 
in agreement with the detection fraction found by Elbaz 
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Figure 2. Detection fraction as a function of Sersic index (upper 
left panel), effective radius (upper right, units of kpc), mass (bot- 
tom left) and redshift (bottom right). The black dots indicate the 
24 fim detection fraction with Poissonian error bars, while the 
red dots show the fraction of 250 fim detected sources. 



et al. (2010) analysing the whole GOODS-N data sample 
with PACS and/or SPIRE detections. Most of the SPIRE 
detected sources in our sample are disk-like objects. More 
precisely, the detection fraction of galaxies with n ^ 2, is 45 
per cent, while we detect only ~ 13 per cent of the spheroid- 
like galaxies (n > 2). Among the MIPS 24 fim detections, 
the detection fractions are in agreement with other literature 
results (e.g. PG08) 

• for both the MIPS- and SPIRE-detected fractions (see 
Figf2J) a trend is present indicating that larger, massive, 
disk-like galaxies are better detected with respect to small, 
spheroidal, less massive galaxies. This would confirm the 
findings of Kriek et al. (2009) who find, using optical spec- 
troscopy, that large disk-like galaxies are much more star- 
forming than compact spheroid-like galaxies 

• as an additional proxy of the compactness of galaxies 
we have split our sample based on galaxy surface density. We 
find that on average less compact galaxies are more easily 
detected with respect to more compact objects. We also find 
a small indication of a larger mean size and mass for detected 
galaxies with respect to non-detections. This is supported by 
the trends shown in Fig. 2 

• by means of fitting to the IR part of the spectral en- 
ergy distribution, we have derived the SFR and sSFR for 
the SPIRE and MIPS sub-samples including the contribu- 
tion from non-detections for SPIRE and accounting for the 
systematic effect of overestimation of the LIR when using 
MIPS24 based measurements. The two estimates result in 
good agreement 

• disk-like galaxies show a mean higher value (~ 3 times) 
of SFR and sSFR with respect to spheroid-like galaxies for 
both SPIRE and MIPS based measurements. Our estimated 
mean sSFR for the two morphological classes at z ~ 2.3 well 
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PACS observations for GOODS North and South fields pro- 
vided by the PACS Evolutionary Project. 
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Figure 3. Comparison of the mean sSFR at z ~ 2.3 (correspond- 
ing to the median redshift of our sample) for SPIRE detected 
galaxies (black square for the global population, red and blue for 
disk- and spheroid-like galaxies, respectively) with other results 
on the evolution of the sSFR. The magenta and cyan shaded re- 
gions show the sSFR measured by PG08 for early and late-type 
galaxies detected at 24/xm, respectively. The red dashed line and 
the blue dotted one are their fitting relations for the redshift evo- 
lution of the sSFR (see text for details) . The light- and dark-green 
lines are the fitting relations for late- and early-type galaxies, re- 
spectively, in Oliver et al. (2010a). The yellow region indicates 
the estimate of Rodighiero et al. (2010), using PACS data. 



match the sSFR evolution of previous studies (PG08, Oliver 
et al. 2010a, Rodighiero et al. 2010). 

The debate on the nature of these different population of 
galaxies (compact and large galaxies) is still open. It is pos- 
sible that they are linked through dissipative major merg- 
ers, as suggested by e.g. Ricciardelli et al. (2010), or they 
might follow two separate evolutionary paths and the dif- 
ferences in the structural parameters simply reflect different 
formation epochs. In the first case, massive galaxies would 
go through a phase of high star formation and the transfor- 
mation from isolated, gas-rich disk-like galaxies with typical 
sizes of ~ 2 — 3 kpc, into compact (R e ~ 1 kpc) less star- 
forming galaxies is driven by violent major merger events, 
compatible with the scenario depicted by theoretical mod- 
els. In the second case, massive galaxies at z > 2 must have 
formed very quickly, and consequently their high stellar den- 
sities could reflect the high gas densities in the primeval 
Universe. The high fraction of merging systems suggested 
by the inspection of ACS- and NICMOS-fiST images for 
the objects presented in this work would go in the direc- 
tion of supporting the first proposed scenario, but further 
investigations are necessary to confirm this hypothesis. 

This work represents the first step in the direction of 
understanding FIR properties of these high- z compact mas- 
sive galaxies and will be followed by a more extensive study 
exploiting the whole GNS sample and including the deep 
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